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ABSTRACT 

We report the discovery of the optical and radio afterglow of GRB 010921, 
the first gamma-ray burst afterglow to be found from a localization by the High 
Energy Transient Explorer (HETE) satellite. We present optical spectroscopy of 
the host galaxy which we find to be a dusty and apparently normal star-forming 
galaxy at z = 0.451. The unusually steep optical spectral slope of the afterglow 
can be explained by heavy extinction, Ay > 0.5 mag, along the line of sight to 
the GRB. Dust with similar Ay for the the host galaxy as a whole appears to be 
required by the measurement of a Balmer decrement in the spectrum of the host 
galaxy. Thanks to the low redshift, continued observations of the afterglow will 
enable the strongest constraints, to date, on the existence of a possible underlying 
supernova. 
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Subject headings: gamma rays: bursts 

1. Introduction 

The High Energy Transient Explorer (HETE-2) was successfully launched and deployed 
on 2000 October 9 (Ricker et al. 2000). As the first satellite entirely dedicated to the detection 
and study of gamma-ray bursts (GRBs), HETE's primary mission is the localization of GRBs 
through their prompt emission and the rapid relaying of the coordinates to ground-based 
observers. Here, we present the discovery of the afterglow of GRB 010921 (Ricker et al. 
2001b), the first afterglow to be identified through localization by HETE-2. 

2. Discovery of an Optical and Radio Transient 

GRB 010921 triggered the FREGATE instrument on board HETE-2 on 2001 September 
21.21934 UT as HETE trigger number 1761 (Ricker et al. 2001a). The GRB was detected 
by only one of HETE's two Wide-field X-ray Monitors, and so localization was only possible 
in a single dimension. Fortunately, the Interplanetary Network was able to provide comple- 
mentary localization and the resultant 310-squarc-arcmin error box was speedily reported 
to the community via the GRB Coordinate Network (GCN) Circulars (Hurley el al. 2001). 
We observed this error box with the Large Format Camera (LFC) on the Hale 200-inch 
telescope at Palomar Observatory, commencing 22 hours after the GRB, in the Sloan Digital 
Sky Survey (SDSS; York et al. 2000) r' filter. Despite the 26-arcmin field-of-view (FOV) of 
the LFC, three pointings were required to cover the entire error box. 

A quick comparison of the images with the Digital Palomar Observatory Sky Survey 
(Djorgovski et al. 1999; also in prep.) failed to identify an afterglow candidate within the 
large error box (Fox et al. 2001). Hoping to identify the afterglow by its variability (i.e., 
decay), we re-observed the error box with LFC five days later. 

Reduction of data from LFC is comphcated by the instrument's large FOV which leads 
to substantial optical distortions in the focal plane, some of which can be removed in software, 
but this does not correct for the inevitable coma at the edges of the field. Furthermore, this 
problem is made worse by our choice of on-chip binning (which however reduced the file 
size of the images). We used the MSCRED package in the Image Reduction and Analysis 
Facility (IRAF)^^ to bias-subtract, flat- field and combine the individual images. 



IRAF is distributed by the National Optical Astronomy Observatories. 
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Comparison of the first- and second-epoch images was done through PSF-matched image 
subtraction (Alard 2000). Subtraction of images was rendered difficult due to the binning of 
pixels. Nonetheless, we identified source A (0:2000 = 22'^55™59**90, ^2000 = 

+40°55 52 9) and 

source B (0:2000 = 

221^56°^ 15*^09, ^2000 = +41°04'48"l), both of which were within the error 
box and were clearly variable. These positions are derived relative to USNO-A2.0 and are 
uncertain by approximately 0.6 arcsec. 

In order to further investigate these two sources, we re-observed the field on 2001 October 
12th and 17th, 21 and 26 days after the GRB respectively. Photometry from these epochs 
revealed that source B had increased in brightness, while source A had faded further and 
appeared to be settling to a constant flux level; see Figures 1 and 2. Source B, which is 
faintly visible at the plate hmit of the second Digitised Sky Survey^^, is likely a variable 
star or some other object unrelated to the GRB, while the light curve of source A bears a 
similarity to the afterglows of other GRBs which are observed to undergo a power-law decay 
followed by leveling off at the flux level of the host galaxy. 

Optical photometry was performed relative to stars in the held calibrated through ob- 
servations of numerous Landolt standards observed at the USNO Flagstaff Station 1.0-meter 
telescope on 2001 September 22. We estimate that these cahbrations have a zero-point error 
of less than 0.03 mag. Reference star magnitudes in the SDSS photometric system were 
derived through application of the transformations of Fukugita et al. (1996). The resulting 
photometry of source A can be found in Table 1. 

A program of observations at the Very Large Array (VLA)^^ were begun commencing 
on 2001 October 17; see Table 2. All radio observations were obtained in the standard 
continuum mode with 2 x 50 MHz contiguous bands. We observed the extra-galactic source 
J2255+420 for phase calibration, and 3C48 (J0137+331) and 3C286 (J1331+305) for flux 
calibration. All data were reduced using the Astronomical Image Processing System (AIPS). 

Coincident with source A, we identifled a radio source at 02000 = 22'^55™59'^931 ±0.018 
and ^2000 = +40°55'52"23 ± 0.20 and with a spectral slope of (3 = 0.35 ± 0.19 (where 
fi, (X u^) between 4.86 and 22.5 GHz. This value of (3 is consistent with that expected from 
the synchrotron emission of an afterglow (/3 = 1/3; Sari, Piran & Narayan 1998), but does 
not exclude the possibihty that the source is a low- luminosity AGN. Further observations 
indicate variability in the source (Table 2), which may be interpreted as due to interstellar 



^^The Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Govern- 
ment grant NAG W-2166. 

^^The National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. NRAO operates the VLA. 
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scintillation, a phenomenon commonly seen in radio afterglows (Frail et al. 1997). 

3. The Afterglow 

As noted above, the optical light curve, a decline followed by leveling off, the variable 
radio emission and the flat radio spectrum, support the identification of source A as the 
afterglow of GRB 010921. However, the data are sparse and the possibility remains that 
source A could be an AGN since AGNs also possess power-law spectra and can be strongly 
variable. However, in section 4 wc note that the optical spectrum of the host docs not show 
any evidence for a central engine. Furthermore, by comparing Palomar LFC images taken on 
2001 September 22 (dominated by the optical transient) with those taken on 2001 October 
17 (dominated by the galaxy), wc find that the transient is offset from the light ccntroid of 
the host galaxy by 0.351 ± 0.049 arcsec, or 2.18 ± 0.30 kpc at the distance of the galaxy 
(see §4). Such angular and physical offsets with respect to their host galaxy arc typical of 
GRBs (Bloom, Kulkarni & Djorgovski 2002). Thus from the totahty of the radio and optical 
data we concluded that source A is not an AGN but rather the afterglow of GRB 010921 
and announced it as such (Price et al. 2001b). 

After accounting for Galactic extinction towards this direction {E(^b-v) = 0.148; Schlegel, 
Finkbeiner & Davis 1998) we fit the optical photometric measurements to a power-law (af- 
terglow) model, Fi, oc We obtain marginally-acceptable fits (x^/DOF = 14.5/10), with 
a = -1.59 ±0.18 and P = -2.22 ±0.23. Park et al. (2002) report similar decay values based 
on observations carried out with a number of telescopes. 

The temporal decay slope is consistent with those seen in previous optical afterglows, 
but the spectral slope is quite steep, suggestive of extinction in the source frame. Adopting 

standard afterglow theory (Sari, Piran & Narayan 1998), wc can deduce the intrinsic spectral 
slope from the temporal data. However, this inference depends on the details of the afterglow 
model (ambient gas, homogeneous versus inhomogeneous; the electron distribution power law 
index, p; the shape of the explosion, spherical versus collimated outflow). The extinction 
along the line-of- sight to the GRB, parametrized by Ay in the rest frame of the source, 
can then be obtained from the deviation of the observed broad-band spectrum from that 
expected from the afterglow model. 

The actual value of the derived Ay is somewhat dependent on the choice of both the ex- 
tinction curve (Milky Way, LMC or SMC), and the afterglow model (collimated or isotropic 
ejecta; cooling break, redward or blueward of the optical band over the course of the obser- 
vations) but, as can be seen from Table 4 it is clear that Ay is in excess of 0.5 mag. The 
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best fit model has a spherical outflow with the cooling frequency below the optical band by 
day 1 and an LMC-hke extinction curve (x^/DOF = 15.5/11). 

4. The Host Galaxy 

Using the Double Spectrograph (Oke & Gunn 1982) at the Cassegrain focus of the 
Hale 200-inch telescope we obtained a spectrum of the host galaxy (Figure 3). We used 
a 1-arcsec wide long-sht, at a position angle 132°, close to the parallactic angle, and a 
158 lines mm~^ grating, giving an effective instrumental resolution FWHM ~ 12A, and an 
approximate wavelength coverage 5.000-10,000 A on the red side. We obtained 2 exposures 
of 1800-s each, with a small dither along the slit. Exposures of arc lamps were obtained for 
calibrations and small wavelength shifts due to the instrument flexure were removed by fitting 
the wavelengths of night sky lines. The net resulting wavelength scatter is ~ 0.4 A. Flux 
cahbration was made through observations of the spectrophotometric standard BD-|-17°4708 
from Oke & Gunn (1983). 

Several emission lines are detected (see Table 3), corresponding to [O II]A3727, H/3, 
[O III]AA4959,5007 and Ha, at a mean redshift of ^ = 0.4509 ± 0.0004. We note that 
the widths of the lines, ranging from 10 A to 13 A, are consistent with the instrumental 
resolution and thus there is no reason to invoke an AGN component to this host galaxy. 
Furthermore, the line flux ratios F5oo7/-^486i a-nd ^37277-^5007 place the host galaxy squarely 
in the range corresponding to normal HIT galaxies on the identification plot of Baldwin, 
Phillips & Terlevich (1981). 

In estimating the properties of the host galaxy, we adopt a standard fiat cosmology with 
Ho — 65 km s~^ Mpc~^, flu — 0.3 and Qa — 0.7. From the observed Balmer decrement and 

an LMC extinction curve, we calculate an optical depth at H/3 due to dust of 1.51 (assuming 
an intrinsic Ha to H/? fiux ratio of 2.87; Mathis 1983); note, we have not accounted for Balmer 
absorption lines in the host galaxy. The corresponding Ay which applies to the galaxy as a 
whole is 1.3 mag. This value of Ay is comparable to the extinction along the line-of-sight to 
the GRB afterglow and thus most likely the extinction observed in the afterglow is due to 
interstellar (as opposed to circumburst) dust. 

There are two other GRB host galaxies for which Balmer line fiuxes have been published, 
the hosts of GRBs 980703 (Djorgovski et al. 1998) and 990712 (Vreeswijk et al. 2001), for 
which we calculate optical depths at H/? of approximately 0.1 and 2.41, respectively, using 
the above method. We note that the Balmer decrement observed for the host of GRB 010921 
falls squarely within the range observed for these other two GRB hosts. 
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We now infer the metallicity (oxygen abundance) using the empirical relation between 
oxygen abundance and log7?23 = log {{[OH] + [0III])/Hf3) (Kobulnicky, Kennicutt & Piza- 
gno 1999). We find logi?23 = 1-05 ± 0.18 which is close to the knee of this relation and thus 
yields an unambiguous measure of the oxygen abundance, 12 + log (0/H) of 8.35 ±0.25. The 
rest-frame equivalent width of the [O II] line, 34.5 ±3.5 A, is typical for field galaxies (Hogg 
et al. 1998), which have equivalent widths spanning - 50 A. Furthermore, using the [O 
II] and Ha line fluxes and the empirical star-formation rate (SFR) estimators of Kennicutt 
(1998) we derive extinction-corrected SFR of 8 Mq yr~^ and 6 Mq yr"^ respectively. All in 
all, we conclude that GRB 010921 occured in a normal, star-forming galaxy with an oxygen 
abundance approximately one-quarter solar. 

We also observed the host galaxy with the Near InfraRed Camera (NIRC; Matthews 
& Soifer 1994) on the Keck I telescope on 2001 November 1 UT in the JHKg bands under 
photometric conditions. The individual frames were flat-fielded and then sky-subtracted and 
combined using the DIMSUM package in IRAF (Eisenhardt et al. 1999). Flux calibration 
was performed through observation of Persson et al. (1998) standards SJ9126 and SJ9134. 
We estimate that the calibration is accurate to within 0.05 mag. The resulting magnitudes 
of the host galaxy, measured within a 1.5-arcsec aperture, are also displayed in table 1. From 
the fiat spectrum in the NIR, we estimate that the stellar population of the host is less than 
1 Gyr in age. 

From the FREGATE 8-400 keV fiuence of 1.5 x 10"^ erg/cm^ (Ricker et al. 2001b), we 
derive the isotropic-equivalent prompt energy release (without any k correction) of 9.0 x 10^^ 
erg. However, it is now generally appreciated that GRBs are beamed or coUimated (Frail 
et al. 2001) and thus the isotropic-equivalent value is an upper limit. Unfortunately, the 
paucity of the afterglow data do not place significant constraint on the opening angle of the 
explosion (see Table 4). 

5. Conclusion 

Here we report the first afterglow resulting from an HETE localization of a GRB. The 
afterglow in itself is similar to afterglows studied to date. However, the GRB occured in a 
low-redshift, z — 0.451, dusty star-forming galaxy. Such low redshift GRBs are valuable for 
two reasons. First and foremost, such low-redshift events offer the best opportunity to study 
the explosion physics including the opportunity to constrain the presence of (or detect) any 
underlying supernova. Second, the low redshift allowed us to determine the metallicity of 
the host galaxy with relative ease. 
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It is perhaps significant to note that the afterglow was discovered through the use 
of image differencing. We suggest that the use of this technique may be a more robust 
method of identifying GRB afterglows than the traditional manual comparison with sky 
survey images, since it enables the detection of an afterglow superimposed on a bright host 
galaxy. This is particularly important for low redshift GRBs. We wonder whether the 
traditional approach — looking for an isolated transient — could be the cause of failure to 
identify optical afterglows in some GRBs. 

We thank Martha Haynes for the generous donation of her Hale 200-inch time without 
which the discovery of this afterglow would not have been possible. We thank Rick Burrus 
for his help with 200-inch observations. PAP gratefully acknowledges an Alex Rodgers 
TraveUing Scholarship. DAF thanks Caltech for the hospitality shown by Caltech during his 
sabbatical leave. SRK and SGD thank NSF for support of our ground-based GRB observing 
program. JSB acknowledges support from the Hertz Foundation in the form of a fellowship. 
KH is grateful for Ulysses support under Gontract JPL 958059, and for HETE support under 
Gontract MIT-SG-R-293291. 
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Table 1. Optical observations of the afterglow of GRB 010921. 



Date (UT) Filter Magnitude Telescope Lead Observer 



2001 


Sep 


22. 


144 


r' 


19 


.593 


± 


0. 


,033 


P200 + LFC 


Haynes/Burrus 


2001 


Sep 


22. 


148 


r' 


19, 


.601 


± 


0. 


,039 


P200 + LFC 


Haynes/Burrus 


2001 


Sep 


22. 


,289 


I 


18, 


.854 


± 


0. 


,064 


USNOFSl.O 


Henden 


2001 


Sep 


22. 


293 


R 


19 


.669 


± 


0. 


,060 


USNOFSl.O 


Henden 


2001 


Sep 


22. 


298 


V 


20, 


.296 


± 


0. 


,070 


USNOFSl.O 


Henden 


2001 


Sep 


22. 


304 


B 


21, 


.159 


± 


0, 


,088 


USNOFSl.O 


Henden 


2001 


Sep 


22. 


313 


U 


20, 


.910 


± 


0, 


,205 


USNOFSl.O 


Henden 


2001 


Sep 


22. 


,321 


R 


19, 


.790 


± 


0, 


,069 


USNOFSl.O 


Henden 


2001 


Sep 


27. 


,354 


^1 


21, 


.620 


± 


0. 


,036 


P200 + LFC 


Djorgovski 


2001 


Oct 


12. 


,226 


R 


21, 


.693 


± 


0. 


,053 


P200 + LFC 


Bloom 


2001 


Oct 


17. 


,145 


r' 


21, 


.994 


± 


0. 


,028 


P200 + LFC 


Mahabal 


2001 


Oct 


17. 


,165 


i' 


21, 


.660 


± 


0. 


,022 


P200 + LFC 


Mahabal 


2001 


Oct 


18. 


,088 


r' 


21, 


.918 


± 


0, 


,030 


P200 + LFC 


Mahabal 


2001 


Oct 


18. 


,110 


i' 


21, 


.670 


± 


0, 


,021 


P200 + LFC 


Mahabal 


2001 


Oct 


19. 


,178 


g' 


22, 


.947 


± 


0, 


,043 


P200 + LFC 


Mahabal 


2001 


Oct 


19. 


,109 




21, 


.977 


± 


0, 


,028 


P200 + LFC 


Mahabal 


2001 


Oct 


19. 


,130 


%' 


21, 


.650 


± 


0. 


,030 


P200 + LFC 


Mahabal 


2001 


Oct 


19. 


,149 


z' 


21, 


.419 


± 


0. 


,043 


P200 + LFC 


Mahabal 


2001 


Oct 


19. 


,253 


B 


23, 


.423 


± 


0. 


,079 


P60 


Fox 


2001 


Oct 


19. 


,206 


V 


22, 


.324 


± 


0. 


,057 


P60 


Fox 


2001 


Oct 


19. 


,272 


R 


21, 


.807 


± 


0. 


,051 


P60 


Fox 


2001 


Nov 


01 


.287 


Ks 


19, 


.069 


± 


0, 


,037 


KI + NIRC 


Kulkarni 


2001 


Nov 


01 


.345 


H 


19, 


.750 


± 


0, 


,035 


Kl + NIRC 


Kulkarni 


2001 


Nov 


01 


.316 


J 


20, 


.338 


± 


0. 


,019 


KI + NIRC 


Kulkarni 



Note. — These measurements have not been corrected for Calactic extinc- 
tion. The error in the measurement includes both statistical and systematic 
errors. A 2.5-arcsec aperture was used for optical measurements, and a 1.5- 
arcsec aperture for the NIR measurements. Telescopes are: P200 — Hale 
Palomar 200-inch; USNOFSl.O — USNO Flagstaff Station 1.0-metre; P60 — 
Palomar 60-inch; KI — Keck 1. 



- 11 - 



Table 2. Radio Observations of GRB 010921 made with the Very Large Array. 



Epoch z/q S±(T 

(UT) (GHz) (/.Jy) 



2001 Oct 17.15 4.86 188 ± 25 

2001 Oct 17.15 8.46 222 ± 16 

2001 Oct 17.15 22.5 330 ± 90 

2001 Oct 18.23 8.46 229 ± 22 

2001 Oct 19.02 4.86 100 ± 28 

2001 Oct 24.02 22.5 37 ± 95 

2001 Oct 26.04 4.86 91 ± 36 

2001 Oct 28.99 4.86 140 ± 28 

2001 Oct 28.99 8.46 158 ± 26 

2001 Oct 29.99 8.46 123 ± 26 



Note. — The columns are 
(left to right), UT date of 
the start of each observation, 
observing frequency, and peak 
flux density at the best fit posi- 
tion of the radio transient, with 
the error given as the root mean 
square noise on the image. 
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Table 3. Lines identified in the spectrum of the host galaxy of GRB 010921. 



Aobs (A) 


Line 


-^obs 


EW (A) 


GW (A) 


TP 

^ corr 


5408.77±0.57 


[0 II] 


2.42±0.23 


50.1±5.1 


11.5±1.1 


24.4±2.3 


5613.6 


[Ne III] 


< 0.98 


< 13.8 




< 8.9 


7051.35±1.00 


H/3 


0.69±0.11 


13.2±2.3 


10.2±1.9 


3.12±0.50 


7192.81±1.76 


[0 III] 


0.41±0.13 


9.2±3.0 


10.8±3.3 


1.75±0.55 


7264.48±0.55 


[0 III] 


2.18±0.14 


43.5±3.5 


12.9±1.1 


9.11±0.59 


9522.46±1.22 


Ha 


3.29±0.59 


52±13 


12.3±3.7 


8.9±1.6 



Note. — Left to right, the columns are the observed wavelength of the 
line, line identification, observed flux corrected for Galactic extinction 
using E(^B-v) = 0.148, observed equivalent width, observed Gaussian 
width and flux corrected for extinction in the host galaxy using an optical 
depth at H/3 of 1.51 and an LMC extinction curve. Fluxes are in units 
of 10~^^ erg/cm^ /s. Given the suggestion that GRB host galaxies may 
exhibit strong [Ne III] emission (Bloom, Djorgovski & Kulkarni 2001) 
we have included an entry for this line. 
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Table 4. Results of extinction curve fits to the optical afterglow photometry. 



Afterglow- 


CoUimated? 


MW (C2 = 


2/3) 


LMC (C2 = 


4/3) 


SMC (c2 = 7/3) 


Model 


yes/no 


Ay (mag) 




Ay (mag) 




Av (mag) 


v < Uc, p = 3.1 


no 


1.05 


18.8 


0.73 


16.1 


0.95 16.1 


V > p = 2.8 


no 


0.76 


16.8 


0.52 


15.5 


0.66 15.5 


V > p — 1.6 


yes 


1.29 


21.0 


0.89 


16.7 


1.17 16.8 


V < Uc-i P — 1-6 


yes 


1.75 


26.0 


1.20 


17.7 


1.60 18.0 



Note. — The extinction curve formulation of Reichart (1999) was used in performing 
these fits. The models allow for both coUimated {tjet < 1 d) and isotropic {tjet > 25 d) 
ejecta, and for the cooling break to be either redward (z/ > or blucward (z/ < z/^) of 
the optical bands over the period of observations of the afterglow (1 d < t < 25 d). Ay 
is source frame V-band extinction; typical errors are 0.15 mag (MW). 0.10 mag (LMC) 
and 0.20 mag (SMC), not including co- variance. Each fit had 11 degrees of freedom. 
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Fig. 1. — Identification of the optical afterglow of GRB 010921. Left and middle: r' images 
from the Hale 200-inch with the LFC on 2001 September 22 and 27 respectively. Right: 
Corresponding image from the DSS-2. Each image is 1.5-arcmin on a side, with North to the 
top and East to the left. The position of the optical afterglow is circled. Star Z (labeled), 
at coordinates 22^56"00'^2 +40°55'22" (J2000) is 2.65" E and 30.35" S of the afterglow, and 
for which we measure R — 21.169 ± 0.043 mag and / — 19.537 ± 0.025 mag. 
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Fig. 2. — The optical light curve of the afterglow of GRB 010921. The displayed fluxes have 
been corrected for Galactic extinction using E(^b-v) — 0.148. The lines indicated the best-fit 
power-law plus galaxy model. 
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Fig. 3. — The spectrum of the host galaxy of GRB 010921. The spectrum has been smoothed 
with a Gaussian of FWHM = 12 A, approximately the instrumental resolution. The contin- 
uum longward of 8500A suffers from poor sky subtraction. Emission lines corresponding to 
[O II], B.p, [O III] and Ho; (labeled) are clearly detected, corresponding to a mean redshift 
of ^ = 0.4509. 



